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Abstract 
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' We present results for two-loop QCD corrections to the short distance coefficient that gov- 

^ ■ erns leptonic decay of Bc-meson. After brief discussion of the technics of asymptotic expansion 

d ', used to perform this calculation, we give comments on the numerical value of two-loop correc- 

I tions and its impact on NRQCD description of i?c-™eson bound state. 
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> '. 1 Introduction 
(N . 

■ 

■ The description of various aspects of Bc-iaeson physics has received recently a lot of attention 
CN i both from theoretical and experimental sides. 5c-meson is the only meson within SM, which 

' contains two heavy quarks of different flavor. For this reason one might expect, that its 

■ spectroscopy, production mechanisms and decays diff'er significantly from those of charmonium, 
, bottomonium as well as hadrons with one heavy quark. Indeed, contrary to the case of J/^ 

Q^l and T-mesons, i?c-meson is a long living system, decaying trough electroweak interactions. 

Q^' To produce 6c-system in fiavor conserving collisions, an additional heavy quark pair should be 

(D . created, what makes the consideration of production mechanism complex already at the leading 

order in the coupling constant. The bound state heavy quarks dynamics, being very similar to 
' that of other quarkonia, at the same time has some differences. With respect to the hierarchy 
^ . of dynamical scales, i?c-meson stands among the families of charmonium cc and bottomonium 

^ I bb and thus could be used to study both quantitatively and conceptually existing effective low 

- - ' energy frameworks for the description of bound state heavy quark dynamics, like NRQCD [2], 

pNRQCD [3] and vNRQCD [4] . So, the study of Bc-meson gives us an opportunity to test SM 
in a completely new dynamical regime, explore it at a different angle, what may give us an 
additional insight into our understanding of SM in general. 

Today, we already have many theoretical predictions, concerning various Bc-uieson prop- 
erties: spectroscopy [5], production [6] and decay channels [7]. In the present work we will 
be interested in radiative corrections to the leptonic decay constant of Bc-meson, entering 
expression for its purely leptonic width. This particular decay mode of Bc-meson is not really 
interesting for a direct measurement — the necessarily spin-fiip makes its branching fraction 
very small. However, the same quantity appears in the formulation of two-point QCD sum 
rules [8]. And the latter may be used to extract an additional information about parameters 
of SM, namely the heavy quark masses. A first step towards systematic formulation of NNLO 
-Bc-meson sum rules, treating Bc-meson as a Coulomb system, was made in [9], where NLO 



^For a review see [1] 
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analysis was presented. In order to describe Sc-meson QCD sum rules at NNLO, one needs 
both 2-loop matching Wilson coefficient of QCD current with Be quantum numbers on cor- 
responding NRQCD current and two point NRQCD correlator, computed at NNLO. Below 
we will present the results of matching procedure. These results are also interesting from the 
point of view of further development of computation technics used in dimensionally regulated 
NRQCD [10]. Here we are dealing with threshold expansion with two mass scales nib and rric, 
what is an extension of previously obtained results for matching of vector current with two 
equal heavy quark masses [11, 12]. The presence of two mass scales of course makes problem 
more complicated, however as we will see the whole calculation could be performed within 
the general framework of threshold expansion in dimensionally regulated NRQCD [10]. Here 
wc would like to stress, that wc treat Sc-nicson as a Coulomb system. So, in our matching 
procedure c-quark is considered to be heavy with respect to factorization scale rric > //fac- 
As a consequence, it is impossible to make connection of our result with the results on the 
matching of heavy-light currents [13]. On the other side the limit rrib = rric is well defined in 
our approach and in the case of vector current matching with equal quark masses we reproduce 
already known results of [11, 12]. 

This paper is organized as follows. In section 2 we introduce necessary notation and 
present one-loop results obtained previously [14]. Section 3 contains two- loop results for Bc- 
mcson leptonic decay constant. Here we also describe shortly various steps of the performed 
calculation. And finally, in section 4 we present our conclusion. 

2 Be - meson leptonic constant: one loop result 

The Bc-meson leptonic decay proceeds through a virtual VF-boson. The VF-boson couples 

to the Be through the axial-vector part of the charged weak current. All QCD effects, both 
perturbative and nonperturbative, enter into the decay rate through the decay constant Jb^, 
defined by the matrix element 

{0\br^5c\BeiP)) = ifB^P^ (1) 

where \Bc{P)) is the state of Sc-uieson with four-momentum P. It has the standard covariant 
normalization {Bc{P')\Bc{P)) = {2'k)'^2E 5^ {P' — P) and its phase has been chosen so that fs^ 
is real and positive. 

To separate short-distance (~ l/mq) and long-distance l/rriQ) QCD contributions to 
fSc we use the formalism of NRQCD [2]. As a result the decay amplitude is factored into 
short-distance coefficients multiplied by NRQCD matrix elements. The Wilson short-distance 
coefficients could be then calculated further order by order in perturbation theory by matching 
a perturbative calculation in the full QCD with the corresponding perturbative calculation in 
NRQCD. It is the subject of this paper to calculate Wilson coefficient in front of NRQCD 
current up to 2-loop order in perturbative QCD. 

Let us now discuss the matching procedure in some detail. We start with NRQCD la- 
grangian describing Bc-meson bound state dynamics 

>Cnrqcd = Aight + i'l {iDo + DV(2me)) Vc + xl {iDo - T>V{2mb)) Xb + ■ ■ ■ , (2) 

where £iight is the usual relativistic lagrangian for gluons and light quarks. The 2-component 
field tpc annihilates charm quarks, while Xb creates bottom antiquarks. The typical relative 
velocity v of heavy quarks inside the meson provides a small parameter that can be used as a 
nonperturbative expansion parameter. 
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To express the decay constant /b^ in terms of NRQCD matrix elements, we express the 
axial-vector current bjf^j^c in terms of NRQCD fields. Only the = component contributes 
to the matrix element (1) in the rest frame of the i?c-meson. This component of the current 
has the following operator expansion in terms of NRQCD fields: 

67°75c = Co{mb,mc) xl'^c + C2{mb,mc) (Dxb)^ • DV'c + (3) 

where Cq and C2 are Wilson short-distance coefHcients that depend on the quark masses rrib 
and nic- By dimensional analysis, the coefficient C2 is proportional to l/rriQ. The contribution 

to the matrix clement {0\b-y^j5c\Bc) from the operator (Dxb)^ ''^''Pc is suppressed by v'^ relative 
to the operator xlV-'o where v is the relative velocity of heavy quarks inside i?c-meson. The 
dots in (3) represent other operators whose contributions are suppressed by higher powers of 

The short-distance coefficient Cq and C2 can be determined by matching perturbative 
calculations of the matrix elements in the full QCD and NRQCD. A convenient choice for 
matching is the matrix element between the vacuum and the state \cb) consisting of a c and a 
b on their perturbative mass shells with nonrelativistic four-momenta p and p' in the center of 
momentum frame: p -|- p' = 0. The matching condition is 



(0|67°75c|c6) 



= Co (0|xiV'c|c6) 

QCD 



+ C2 (0|(Dx6)^-DV'c|c6) 

NRQCD 



+..., (4) 

NRQCD 




where QCD and NRQCD represent perturbative QCD and perturbative NRQCD, respectively. 
At leading order in Ug, the matrix element on the left hand side of (4) is p)7''75Uc(p). 
The Dirac spinors are 

ITT' ™ / C \ 

(5) 
(6) 

where and rj arc 2-component spinors and Eq = rriq+'p^. Making a nonrelativistic expansion 
of the spinors to second order in p/mg, we find 

^.(-P)7°75«c(p) « vUc (l - I f!I!±±^)\^ + ..]. (7) 
° \ 8 V mbiric ) ) 

At leading order in a^, the matrix elements on the right hand side of (4) are ryj^c and p^?7^Cc- 
The short distance coefficients are therefore Co = 1 and 

where mred = mbrric/ [rub + rric) is the reduced mass of the system. 

In order to determine the short distance coefficients up to order a^, we must calculate the 
matrix elements on both sides of (4) up to a^. Up to this order it is sufficient to calculate 
radiative corrections only for the coefficient Cq, since the contribution proportional to C2 is 
suppressed by additional power of v^. The coefficient Cq can be isolated by taking the limit 
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p — 0, in which case the matrix element of (Dxh)^ • D^'c vanishes. Such calculation at the 
one-loop order was performed in [14], where the following result was obtained for Cq: 



Co = i + 



IT 



] log — 

rub + rUc rUc 



(9) 



mred is the scale of the running coupling constant. To the accuracy of this calculation, any 
scale of order rrib or nic would be equally correct. 

Thus, the final result for the decay constant of the Be is 

ifB^MB, = Co {0\xIA\Bc) + C2 (0|(DX6)^ -D^elSc). (10) 

with short-distance coefficient Co up to next-to- leading order in as given in (9), while C2 is 
given to leading order in (8). The uncertainties consist of the perturbative errors in the short 
distance coefficients and an error of relative order f ^ from the neglected matrix elements that 
are higher order in v"^. The matrix element (0|x^Vc|-Bc) can be estimated using wavefunctions 
at the origin from nonrelativistic potential models: 

mxl^c\Bc)\' ^ 2MB^^\Rm'- (11) 
The factor of 2Mb^ takes into account the relativistic normalization of the state \Bc). 



3 Two-loop result 

In this section we consider the two- loop contribution to i?c-meson leptonic constant. An 
analogous calculation, but for the case of vector current and equal quark masses was done 
previously in [11, 12]. Here we are following the approach of [11, 12], extending it, where it 
is necessary, to the case of two different mass scales involved in our problem. In the case of 
Sc-meson we consider the matching of QCD axial current on the NRQCD one. 

To compute Wilson short distance coefficient Co up to ct^ order, we replace i?c-meson state 
by a free quark-antiquark pair of on-shell heavy quarks with small relative velocity. In terms of 
the on-shell matrix elements for heavy quarks with small relative velocity: p = qnii / {nii + 1712) 
and p' = qm2/{mi + m2) {q is the momentum of the gauge VF-boson) ^. The matching equation 
has the following form 

Z2,QCD^QCD = CoZ2,NRQCD-^J^NRQCDrNRQCD + ©('W^), (12) 

where Tqcd and rNRQCD are amputated bare axial-vector vertices in QCD and NRQCD, 
correspondingly. Here, we would like to note, that since terms of order v'^ are not needed 
to determine Co (see Eq. (4)), one can set relative momentum to zero and compute Green 
functions directly at the threshold. Note, that here we deal with partially conserved axial- 
vector current and its QCD renormalization constant is given by Zj^qqb = 1, as in the case of 
vector current. 

In order to obtain the coefficient Cq at two loops one has to perform the following steps: 

1) sTim over all contributions to F's, including the one-loop term, multiplied by the two- 
loop QCD on-shell wave function renormalization constant [15], -Z^jqcd ^^'^ one-loop 
NRQCD-current renormalization constant (at one-loop it is equal to unity), 

^This and previous definitions of external heavy quark momenta lead to the same value of the matching coefficient, 
but latter is more convenient in two- loop calculations. 
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Figure 1: QCD Feynman diagrams relevant for the calculation of the two-loop corrections to the 
Be structure constant. 



2) perform the one-loop renormalization of the coupUng and masses. 

After these manipulations the remaining poles in e belong to the anomalous dimension of the 
corresponding NRQCD current and finite part correspond to short-distance constant Cq. 

In step 1) calculating two- loop contribution to vector- axial vertex, it is convenient to write 
it in the following form 



mi + m2 



q^ilbFiiq^) c{p). 



(13) 



The formfactors Fi{q'^) and F2{q^) can be singled out by imposing appropriate projections on 
them and thus our calculation will reduce to the calculation of diagrams with nonvanishing 
numerators. 

To simplify matching calculation at 2-loop order we use threshold expansion in relative 
velocity v of [10]. The latter is obtained by writing down contributions corresponding to hard 
{k ~ mg), soft {k ~ rnqv), potential {ko ~ mQV^,ki ~ rngv) and ultrasoft {k ~ rnqv^) 
regions^. The contributions from the last three regions can be easily identified with NRQCD 
diagrams, that appear in the calculation of rNRQCD- Thus, they simply drop out from the 
matching relation (12) and it will suffice to compute only contribution to threshold expansion 
of Fqcd of the region, where all momenta are hard. Diagrams, whose hard contribution need 
to be computed, are depicted in Fig. 1. Both form factors, introduced above have Coulomb 
singularities at the threshold and diverge as l/v^. However, these singularities appear only in 
the soft, potential and ultrasoft contributions, while the hard contribution is well defined in 
dimensional regularization directly at the threshold. Once the relative velocity is set to zero, the 
corresponding loop integrals are considerably simplified and depend only on two scales M and 
m, where M is the mass of the heaviest quark and ra is the mass of the lighter quark. All loop 
integrals with irreducible numerators can then be further reduced to the scalar integrals with 
increased space-time dimension [16]. For the scalar integrals we derived recurrence relations 
using integration by parts method of [17], which allowed us to reduce all scalar integrals to 
sunsets. The sunset integrals obtained were further reduced to master integrals with the use 



^Here k is the integration momentum. 



5 



of recurrence relations of [18, 19, 20, 21]. We would like to note, that our reduction procedure 
is valid for arbitrary values of heavy quark masses, while to calculate some master integrals 
an expansion in parameter r = m/M was performed. All master integrals are of the sunset 
type discussed in the paper [21]. In particular there we discuss in full details the procedure of 
expansion in parameter r. Here, we describe only shortly the calculation procedure and more 
detailed exposition will be given elsewhere. The results, obtained for each of diagrams in Fig. 
1 could be found in Appendix A. 

After performing the steps of renormalization procedure described above, we have the 
following MS expression for the anomalous dimension of NRQCR current, which first arises 
here at the two-loop order (r = m/M) 



7J,nrqcd 



dlnZj^NRQCD 

dln/Lt 



= - C 



fl 



1+r 



+ Cpc\'^(^)\o{al). (14) 



This scale dependence is compensated by the scale dependence of the two-loop matching co- 
efficient C2{M/ijl), which is defined as 



CoK, -) = 1 + ci(M//.)^^ + C2{M/^^) f ^ 



(M)y 



+ ..., 



(15) 



In order to obtain C2 in analytical form, we perform the expansion of master integrals, con- 
sidering the parameter r = m/M as small (in the practice r = mc/mh ~ 0.3). Below we give 
the result of expansion up to O(r^) inclusive. Separating the different colour group factors, in 
C2(M//x) in the MS we have: 

C2{M/n) = CpC2,A + CfCaC2,NA + CFTFniC2,L + CfTfC2,H, 



C2,A 




^^ + -^ + ^^ + -^ln(2)-^ 
3 32 ^ ^ 48 



+ 



29 
16 



165 



257r2 



3.U + _L+12C3 ^2 



C2,NA 



H + 

16 



ISvr 



127 



4 V/i^ , 

25,2 65, 

+r L-\ L 

16 24 48 



L + L 

24 96 



K 1 2, 4l7r^ 23 

-Cs - -vr^ ln(2) H 

2 ^ ^ 96 48 



27r2 
3 



49 
32 



15 . 23, 197r2 

—L^ L + 

8 16 



31 

48 ^36 



C2,L 



1 



4 

C2,H = T- 



11 , 57r2 1 
24^+144 + 12+^ 



"L2 


l' 




' L 


L2" 


2 


12 


12 


2 



IIL 



24 

L2 



197r2 1147 



96 



288 



L2 
2 



L 
12 



237r^ 
96 



103 
144 



_ 89L 237r^ 11527 
2" + "60"+ 96 + 3600 



Here ^3 = 1-202. . ., M = mfe, r = mc/mi,, L = ln(r) and we have taken all fermions with 
masses less than and mc as massless. The heavy quark masses, appearing in our formulae 
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are pole quaxk masses. All calculations were performed with the help of computer algebra 
system FORM [22]. As a check for our calculations, we employed a general covariant gauge. 
To test our FORM programs we also have reproduced already known results for the sort 
distance corrections to the vector current with equal masses. 

Now let us discuss the expression for 2-loop anomalous dimension of Sc-meson NRQCD 
current. It is an exact formula in a sense, that we resumed a whole series in r. In the limit of 
equal quark masses it differs from the result obtained in the case of vector current. This fact 
is explained by the contribution, coming from Breit-Fermi potential, namely spin-spin and 
spin-orbital heavy quark interaction, which in the case of Sc-meson has a different expression 
compared to that computed for J/* and T-mesons [23]. This result for anomalous expression 
could be also obtained considering NNLO corrections to the Coulomb Green function, given 
by Breit-Fermi and non-abelian potentials. 

Numerically, the scale dependence of the two-loop matching coefficient (16) is large. We 
expect it to be compensated by the scale dependence of 5c-meson wave function computed 
within the lattice NRQCD framework with the scale dependence of the results introduced by 
imposing ultraviolet cutoff/factorization scale. The results of nonrelativistic potential models 
should be considered as being obtained at a scale typical for Bc-raeson bound state dynamics 
~ 1 GeV. For numerical estimates we consider /U = 1 GeV as an adequate bound state scale 
and use its value for factorization scale //fac- In the case of Sc-meson and pole quark masses 
{rrib = 4.8 GeV, rric = 1.65 GeV) we have 



fE 



1 - 1.48 



a. 



- 24.24 



TT 



a. 



TT 



fNR 



(21) 



Here we see that 2-loop corrections are large and constitute pa 100% of 1-loop correction. Even 

less favorable situation was detected earlier for J/'f states [11]. The variation of factorization 
scale does not help to solve this problem. So, at this point one may conclude that application 
of NRQCD to Sc-meson as well as for J/* is not well justified. However, while in the case of 
J/*-meson 2-loop corrections far exceed 1-loop, one may speculate that situation in 5c-meson 
is marginal, in a sense that one may still obtain sensible results using NRQCD framework 
as well as treating Sc-meson as a heavy-light system. This way one may study the change 
in heavy quarks dynamics between bottomonium and charmonium families. Another possible 
reason for such large correction is our use of particular renormalization scheme MS, which 
may turn out not to work well for heavy quark bound state systems. To prove or disprove this 
statement one has to calculate 2-loop corrections to any relation between physical observables, 
from which quarkonium wave function at the origin is eliminated [11]. 



4 Conclusion 

In this paper we presented the results for the two-loop radiative corrections to the short distance 
coefficient in matching between QCD and NRQCD for the i?c-™eson currents. It turns out 
that the numerical value of these corrections are large. There are several conclusions, which 
could be made here. First NRQCD description is inadequate not only for charmonium family 
but also for Sc-meson. To see whether it is true, we suppose to perform a systematical study 
of NRQCD B c-meson sum rules at NNLO in one of our next publications. Second, is it the 
use of MS renormalization scheme, which makes things look nasty?. To answer this question 
one needs to compute some other decay mode of either J/* or Be -mesons at 2-loop level and 
analyze perturbative corrections to the ratio of these decay modes. 
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A 2 loop bare hard contribution 

Here we collected results for hard contributions of diagrams, shown in Fig. 1. Below, to 
save space, we present only coefficient of (a^/vr)^ (e"'^M'^/(47r//^))~^^ for each of the digrams, 
presented above. Though our calculations were performed in general covariant gauge, we 
present here for brevity only the results in the Feynman gauge. 

• diagram 1 (Fig. la) 



9 1 
32i2 +i 



3 

64 



9L' 



11 



4 12 

19 



6 





9L 


,9L 
+ r[- 




4 


~ T 


2 ^ 


L 


-3C3- 


2997r2 


41 


2 ~ 


192 


~ 128 


L + 


491L 
48 


+ 12C3 - 


497r2 
24 



2r_2 



9L, 



o 3 . 9L ^ TT^ 15 
-3L2 + -TT^L - — - 6C3 + — - — 



401 
'36 



}. 



32 

(22) 



• diagram 2 (Fig. Id, with selfenergy for quark with mass M) 

3, 3L 3 ,3L 



^2 f 3 2, 1 /1.3L 3, 3L 3 

32?(i-;:) + 7h:[T + 8l-T-64-'-l 



3t o 9-^ 3- 
+ 4]+^[T + 4- 



1 

r 

+r 



3,3 7r2 13 

-L^ + -LH \ 

8 4 32 12 



9 



13 



+ -L' + —L + 



457r^ 



64 



+ 



IT 



17, 



63 
32 



2041 
384 

2 1001 



5,2 31 , TT 

-L^ H L H \ 

8 6 4 288 



(23) 



• diagram 3 (Fig. Id, with selfenergy for quark with mass m) 



+-L^ + 3L + 



3. 51 ,3, 

--L h r\-L + 

8 64 ^2 



.15 



-L + 



2477r^ 
192 



+ 



1057 



384 



„ ,2 29 , llTT^ 37 

SL^ -\ L -\ \ 

16 96 48 



l^L^ - —L - — - — 
T ~ ~ ~8~ ~ 32 



}. 



(24) 



• diagram 4 (Fig. le, with two gluons coupled to quark with mass M) 



15 1/19 9 9 
64i^ + i 128~ 16^ + '"^16^ 



9 r2 3 ^ 

H L 

16 16 



Ytt^ 329 

\ \-r 

128 256 



16^ +^'^16-16''^; 

15,2 7r2, 11, 3^ 

H L H L C3 

16 6 16 4^^ 



257r^ 
48 



21 
64 



8 



13 



TT 



115 



—L' - —L - —L + -Xz + — - — 



16 



96 



16 144 



}. 



(25) 



diagram 5 (Fig. le, witii two gluons coupled to quark with mass m) 



15 1 
64^2 +i 

15 



9 



_2 
TT 



+— + —L 
8 4 



15^ vr^ 19 
~16 ~ 16 128 ■ ■ '16" ' 16 

2 5^3. 22l7r2 473 



9 vr^ 



16 4^"^ 



H h r 

384 256 



2 



vr^ ^ 3 , 3 , vr^ 21 

— L H L + -Ca H 

4 16 4^^ 12 64 



7r% 83, 3^ vr^ 353 

2L^ H L L - -Cs H \ 

4 96 4^ 8 288 



(26) 



diagram 6 (Fig. Ic, with two gluons coupled to quark with mass M) 
<,HC.-2C.){-i.i(A._- 3 . 3,.^ 3., 5 



16' 



64£2 £Vl6 128 16 16 / 16 " 



16 



vr^ , , 337r2 287 
-T^^('^ + W + 256 +'' 



16^ -^-8^^ + T^^(2)+12-64 



1 

16^ 



11 
48 



57r^ 



71 



— ln(2) + + 

4 ^ 48 36 



(27) 



• diagram 7 (Fig. Ic, with two gluons coupled to quark with mass m) 



rL H HL + -L Cs 

128 16 16 / 8 4 16^ 



5 2, . ^ 777r2 215 
+-7r^ ln(2) + + — + r 



128 256 
"16^ "4^+8^^ 



K2 25, 

-L^ H L 

\ 32 



3 . 3 2- Qtt^ 



-C3-^vr^ln(2) 



53 

32" ~ 32 



VM2) + 



llTT^ 137 

24 "eT 



(28) 



diagram 8 (Fig. lb) 
Cf{Ca-2Cf){ \ 



!L 1 + -TT^L - -L - -Cs - -7rMn(2) - 

8 8^4 ^ ^ 16 



3 2- "-^ 37r 



16 16 y 24 

" 3 n , 5 , 3^ ^ 7r2 59 

— H L L + -Cs H ln(2) \ 

16 24 32 4^"^ 2 ^ ' 32 64 

' 11 ,2 , 11 3, vr^ , , 357r2 

L + — L - -Cs ln(2 

16 12 6 4^ 2 ^ ^ 48 



2 35 
12 



2119 
576 



}, (29) 



• diagram 9 (Fig. If) 
19 



+ 



1/33 19, 19 , 19 2. \ 19 ,2 7 , 957r^ 1315 
- L + —rL - —r^L + —L^ - —L + — — + 



128^2 e V256 32 



32 



32 



16 



32 



768 



512 



— r 



— L2 h L 

16 64 



19,2 5 , 
— L2 h L 

16 64 



(30) 
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diagram 10 (Fig. Ig) 



CfCa\ ^ + L H rL r^L H L'' L H ^ 

128£2 e I 256 32 32 32 / 16 32 768 512 



+r 



1 



1 



—L - —U 



64 



16 



16 



64 



}. 



(31) 



diagram 11 (Fig. Ih, with light quark loop) 



8£2 e I 2 16 2 ^ ' I 48 32 ^ ' 



+ -L 



(32) 



diagram 12 (Fig. Ih, with quark of mass M in the loop) 



^ 4e2 ^ £ 12 48 2 ^ 



+r 



-L' + -L + 



IT 



211 
144 



l-r)L 



-L 



2 13 



L 

3 ~ T ~ 32 
^ TT^ 6007 
30 ~ "e" 3600 



}, (33) 



diagram 13 (Fig. Ih, with quark of mass m in the loop) 

^ ^ r 1 1/, 5 1 x\ „.2 2^ 137r^ 563 



+r 



t2 Ir 77r2 3 

+ -L h - 

4 32 4 



137r2 563 
"9 
137r^ " 
32 



(34) 



where L = Inr and ni is the number of light fermions. 
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